Effective Coulomb Logarithm for One Component Plasma 
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An expression for the efTective Coulomb logarithm in one-component-plasma is proposed, which 
allows to extend the applicability of the classical formula for the self-diffusion coefficient to the 
strongly coupled regime. The proposed analytical approximation demonstrates reasonable agree- 
ment with previous numerical simulation results. Relevance to weakly screened Yukawa systems 
(and, in particular, complex plasmas) is discussed. 
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The one-component-plasma (OCP) model is an ide- 
alized system of point charges (e.g., ions) immersed in 
a neutralizing uniform background of opposite charges 
(e.g., electrons). This model is of considerable interest 
from the fundamental point of view and is relevant to 
a wide class of physical systems, including for example 
laboratory and space plasmas, planetary interiors, white 
dwarfs, liquid metals, and electrolytes. There are also 
relations to various soft matter systems such as charged 
colloidal suspensions and complex (dusty) plasmas. It is 
conventional to characterize the OCP system with the 
ion density n and the temperature T by the dimension- 
less coupling parameter F — e'^ /aT, where e is the charge 
and a = (47rn/3)~^/^ is the Wigner-Scitz radius. In the 
limit of weak coupling, F <C 1, ions form a disordered 
gas-like state. Correlations increase with coupling and, 
at F > 1, OCP exhibits properties characteristic of a 
fluid-like phase. The fluid-solid phase transition occurs 
at Fm ^ 172J. 

The familiar expression for the relaxation frequency in 
ion-ion collisions is^ 



i.0 = (37r)-V2^pr3/2A, 



(1) 



where Wp = ^J^ne^n/m is the plasma frequency, m is the 
ion mass, and A is the so-called Coulomb logarithm. The 
Coulomb logarithm is well defined in the weakly coupled 
regime (F ^ 1). Within the framework of kinetic theory, 
the dominant logarithmic term in A has the form 



A~-ln[l + (ADfc„ 



(2) 



where fc^ax is the upper limit of integration over k 
(corresponding to the exclusion of a close proximity of 
the ion, where kinetic theory is not applicable) and 
Ad = \JT /A-Ke^n is the Debye radius. The proper choice 
of fcmax in this regime is the inverse Coulomb radius 
(_Rc = e^/T). The ion- ion repulsion is very strong for 
r < Rq so that no other ions are essentially present 



in this region and the kinetic approach is inappropriate. 
This results in the well known expression 



A ~ In (AD/i?c) == - ln(%/3F 



3/2\ 



(3) 



The contribution to the relaxation rate from the close 
proximity of the ion can be estimated as vi ^ navr ^ 
n((?lTf^Tlfn ~ WpF^/^ where a ~ Tri?^ ig the rele- 
vant scale of the momentum transfer cross section and 
vt ~ yT/m is the thermal velocity. Since vq/vi ^ A, 
the kinetic approach delivers sufficient accuracy when 
A ^ 1 or, equivalently, F ^ 1. 

Elementary theory of diffusion in weakly coupled 
regime (binary collisions dominate) relates the self- 
diffusion coefficient D to the relaxation frequency v via 
D = T/{mv). Using Eq. ([T|) and conventional normal- 
ization, Dif = D/uj^a?', we obtain 



D^ 



1 



3r5/2A' 



(4) 
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The equations ([3]) and (|4]), which are often referred to 
as the Chapman-Spitzer (CS) result, are essentially ex- 
act in the weakly coupled OCP. More precisely, they pro- 
vide "logarithmic accuracy" in the sense that other terms 
0(1) have been neglected in comparison with the dom- 
inant (large) logarithmic terms in A, as just discussed 
above. The question should naturally arise, whether 
the effective Coulomb logarithm ^ can be meaningfully 
modified in such a way that the applicability of Eq. ^ 
is extended into the regime of strong coupling. 

In this Brief Communication, it is proposed to use 
the approach of Nordholnif^ who showed a simple way 
to greatly reduce inaccuracy of the traditional Debye- 
Hiickel theory with respect to evaluating thermodynam- 
ics properties of moderately and strongly coupled OCP. 
He essentially recognized that the exponential ion density 
must be truncated close to a test ion so as not to become 
negative upon linearization. This resulted in a modifica- 
tion that he called Debye-Hiickel plus hole (DHH) the- 
ory. In the context of the problem addressed here, such 
a hole (void) around a test ion which contains no other 
ions, corresponds to the region where kinetic approach 
should not be applied. Thus, regardless of other possible 
limitations^ a reasonable choice of /cmax would be sim- 
ply the inverse hole radius h~^ . In the weakly coupled 
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FIG. 1: (color online). Reduced self-diffusion coefficient D, 
of the one-component-plasma versus the coupling parameter 
r. Circles: Molecular dynamics results.— Curve: Analytical 
model of Eqs. ©, (O, and 0. 



regime the result remains unchanged, since h ~ Rq in 
this case. In the opposite, strongly coupled regime, we 
have h —^ a, the upper possible limit of the hole radius. 
Using Eq. ([5]) we immediately obtain 



A ~ i In (1 
2 ^ 



i/3r). 



(5) 



In the limit T > 1 this yields A ~ l/(6r) and vq ~ 
(l/6v37r)a;pvr. The contribution to the momentum 
transfer from the excluded region (hole) can be roughly 
estimated considering ion scattering from a hard sphere 



of radius ft ~ a. We get i/i 



navT 



'VT/r 



LUp/y/T. Thus, in the strongly coupled regime, we have 
vq/vi ^ r ^ 1, i.e. kinetic contribution from distant 
collisions should again provide dominant contribution to 
the relaxation rate. 

If this (to some extent naive) picture is relevant, there 
is a hope that a proper interpolation between the weakly 
and strongly coupling limits would adequately repro- 
duce the effective momentum transfer rate ([T|) and self- 
diffusion coefficient (|3]) in a wide range of F. We proceed 
by using the original result by Nordholro^ for the hole 
radius 



ft/Ao = 1 + (SF) 



3/2 



1/3 



1, 



(6) 



and postulating the following expression for the effective 
Coulomb logarithm 



Aeff = \ In [1 + (Ad//i)'] . 



(7) 



Comparison between the present analytical result of 
Eqs. (|4|), (J6|), (l7|) and previous molecular dynamics sim- 
ulation is shown in Figure [TJ The agreement is quite 
reasonable, especially taking into account the approxi- 
mate character of the present model. Although, more 
accurate fits of the self-diffusion coefficient in Yukawa (in- 
cluding limiting case of the OCP) systems have recently 



become available?^ they treat separately weakly coupled 
and strongly coupled regimes, while the present model 
provides reasonable fit over the entire region of coupling, 
up to the fluid-solid transition. 

Moreover, one may expect that the proposed form 
of the effective Coulomb logarithm (possibly with some 
modifications) can be useful in a more general context 
of strongly coupled plasmas. As a demonstration, we 
refer to the results from molecular dynamics simula- 
tion of electron-ion temperature relaxation in a classical 
Coulomb plasma in the regime 0.03 ^ T < 3, described in 
Ref. 0. In this numerical study, positively charged elec- 
trons and ions on a neutralizing (negative) background 
were used to avoid recombination, reassembling classical 
OCP to some extent. A simple numerical fit to the sim- 
ulations results in the entire regime investigated has a 
formi 



A ~ In I 1 + 0.4F"3/2 



(8) 



A comparison between Aeff from Eq. ([7]) and A from 
Eq. ([5]) is shown in Fig. [31 Present approximation ([7]) 
somewhat overestimates the fit ([5]), but the relative de- 
viation (which increases with F) does not exceed ^ 30% 
for the highest F investigated. 

As a final remark, let us briefly discuss the rele- 
vance of the obtained results in the context of complex 
plasmas. Complex (dusty) plasmas consist of highly 
charged micron-size particles in a neutralizing plasma 
background. Signiflcant interest to these systems is in 
large part related to exciting possibilities to use them 
(complimentary to other soft matter systems) to study 
a broad range of important fundamental processes (e.g., 
phase transitions, self-organizations, transport, rheology, 
linear and non-linear waves, etc) as well as with nu- 
merous applications (industrial plasma processes, fusion, 
etc.).— ^— The particles in complex plasmas interact with 
each other via the screened Coulomb (Dcbyc-Hiickcl or 
Yukawa) potential (as long as interparticlc distance does 
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FIG. 2: (color online). Coulomb logarithm A versus the cou- 
pling parameter F. The solid (bluelcurve is a numerical fit to 
the MD simulation data from Ref. [T]. The dashed (red) curve 
corresponds to the effective Coulomb logarithm proposed in 
the present work [Eqs. ©-(0]. 
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FIG. 3: (color online) Reduced self-diffusion coefficient of 
one-component Yukawa systems as a function of the relative 
coupling strength F/Fm- Symbols: Molecular dynamics re- 
sults.— i^ The solid curve corresponds to the approximation 
((9]) for the strongly coupled regime. 



not exceed several screening lengths) J^ In addition, the 
particles are coupled to the surrounding medium (mainly 
to the neutral gas, since the plasma ionization fraction 
is normally quite low) and this coupling can be varied 
in a broad range. At first glance, there can be very lit- 
tle in common between the particle diffusion in complex 
plasmas and that in the OCP. 

Nevertheless, important similarities do exist. To 
proceed further it is convenient to introduce another 
normalization for the self-diffusion coefficient D^ = 
Dn^^^y^m/T, which is often referred to as Rosenfeld's 
normalization^^ It was shown recently, that in weakly 
screened Yukawa systems Dr is a universal function of 
the relative coupling strength T/Fm in the strongly cou- 
pling regimeJ^ Therefore, the proposed approximation 
for the OCP, would be sufficient to describe the diffusion 
in such systems. Using the asymptote A ~ l/(6r), we 
easily get D^ ~ 6.60/F. With Fm — 172 we obtain 



Dn ^ 0.0384 (Fm/F) . 



(9) 



In this form, Eq. (jH]) is applicable not only to the OCP, 
but also to weakly screened Yukawa systems. This is 
illustrated in Fig. [3] The agreement between the sim- 
ple expression ^ and numerical data is quite good in 
the regime F > 20 and k < 2, where k = o/Ad is the 



screening parameter J^ Since quite accurate fits for the 
dependence Fm on k are available^iiii"— Eq. ^ repre- 
sents a useful tool for estimating self-diffusion in weakly 
screened Yukawa systems (e.g. complex plasmas in the 
regime of weak damping) . 

The effect of damping due to the frictional coupling 
between the particles and the neutral gas should also be 
discussed. A simple heuristic expression relating the un- 
damped (Newtonian) and overdamped (Brownian) dif- 
fusivities has been recently proposed;^ In the strongly 
coupled regime, both reduced diffusivities are propor- 
tional to each other, and the relation is particularly sim- 
ple D/Dq ~ ciDr, where D is the actual diffusion co- 
efficient. Do — T / {mvi-[) is the bare Brownian diffusion 
coefficient (no interparticle interactions), v^ is the macro- 
scopic friction rate due to collisions between the particles 
and neutral atoms, and ci = 3.3176 is a numerical con- 
stant i^ This relation can be used to evaluate the particle 
diffusion coefficient in complex plasmas in the limit of 
strong damping (e.g., high neutral gas pressures). More- 
over, it has been suggested that in the regime of interme- 
diate damping, the reduced diffusion coefficient can be 
universally expressed via the dimensionless " damping in- 
dex" ji^ This quantity is defined as the ratio between the 
typical interparticle spacing and mean ballistic free path 
of the particles S, — i^h-/{n-^^^VT) and is a natural mea- 
sure of the damping strength. The expression proposed 
in Ref. [la reads 



D/Da ~ Dn 






(10) 



It is expected to describe particle diffusion in weakly 
screened strongly coupled complex plasmas in the entire 
range of frictional dissipation. 

To conclude, a heuristic expression for the effective 
Coulomb logarithm for the one-component-plasma model 
has been proposed. This expression can be used to de- 
scribe the self-diffusion in the OCP fluid in the entire 
regime of coupling strengths below crystallization. It 
yields essentially exact results in the limit of weak cou- 
pling, renders quite accurate results in the limit of strong 
coupling (near the crystallization point) and provides 
reasonable approximation in the intermediate regime. In 
addition, the obtained results can be of certain value 
in the context of strongly coupled plasmas and weakly 
screened Yukawa systems (complex plasmas). 



S. Hamaguchi, R. T. Farouki, and D. H. E. Dubin, Phys. 
Rev. E 56, 4671 (1997). 

J. D. Huba, NRL Plasma Formulary (Naval Research Lab- 
oratory, Washington, 2004). 

S. Nordholm, Chem. Phys. Lett. 105, 302 (1984). 
These can be for instance associated with the requitements 
to linearize the kinetic equation. This is not discussed in 



detail here, since the plausibility of the final result will 
be checked against "exact" results from numerical simula- 
tions. 

J. P. Hansen, I. R. McDonald, and E. L. Pollock, Phys. 
Rev. A 11, 1025 (1975). 

J. Daligault, Phys. Rev. E 86, 047401 (2012). 
G. Dimonte and J. Daligault, Phys. Rev. Lett. 101, 135001 



(2008). 
^ H. Ohta and S. Hamaguchi, Phys. Plasmas 7, 4506 (2000). 
® M. Chaudhuri, A. V. Ivlev, S. A. Khrapak, H. M. Thomas, 

G. E. MorfiU, Soft Matter 7, 1287 (2011). 
1° V. E. Fortov and G. E. MorfiU (Editors), Complex and 

dusty plasmas: From, Laboratory to Space, (CRC Press 

2010). 
^^ M. Bonitz, C. Henning, and D. Block, Rep. Progr. Phys. 

73, 066501 (2010). 
^2 V. E. Fortov, A. V. Ivlev, S. A. Khrapak, A. G. Khrapak, 

and G. E. MorfiU, Phys. Rep. 421, 1 (2005). 
" S. Khrapak and G. MorfiU, Contrib. Plasma Phys. 49, 148 

(2009). 
" Ya. Rosenfeld, Phys. Rev. A 15, 2545 (1977). 
^^ S. A. Khrapak, O. S. Vaulina, and G. E. MorfiU, Phys. 



Plasmas 19, 034503 (2012). 

^® In complex plasmas, the interparticle separation is often 
characterized by A = n~^'^ and the screening parameter 
is defined as k = A/ Ad. In this notation, the condition 
of weak screening becomes k < 3.2 (A ~ 1.6a), which is 
satisfied under many realistic experimental conditions. 

^■^ O. S. Vaulina and S. A. Khrapak, JETP 92, 228 (2001). 

1* O. Vaulina, S. Khrapak, and G. MorfiU, Phys. Rev. E 66, 
016404 (2002). 

^^ S. A. Khrapak and G. E. MorfiU, Phys. Rev. Lett. 103, 
255003 (2009). 

^" M. J. Pond, J. R. Errington, and M. Truskett, Soft Matter 
7, 9859 (2011). 



